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Abstract 

We present a dynamical model that successfully explains the observed time evolution of the mag- 
netization in diluted magnetic semiconductor quantum wells after weak laser excitation. Based on 
the pseudo-fermion formalism and a second order many-particle expansion of the exact p — d ex- 
change interaction, our approach goes beyond the usual mean-field approximation. It includes both 
the sub-picosecond demagnetization dynamics and the slower relaxation processes which restore 
the initial ferromagnetic order in a nanosecond time scale. In agreement with experimental results, 
our numerical simulations show that, depending on the value of the initial lattice temperature, a 
subsequent enhancement of the total magnetization may be observed within a time scale of few 
hundreds of picoseconds. 

PACS numbers: 
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I. INTRODUCTION 



Ultrafast light-induced magnetization dynamics in ferromagnetic films and in Diluted 
Magnetic Semiconductors (DMS) is today a very active area of research. From the observa- 
tion of the ultrafast dynamics of the spin magnetization in nickel films l| and the analogous 
processes in ferromagnetic semiconductors 2J, special interest has been devoted to the de- 
velopment of dynamical models able to mimic the time evolution of the magnetization on 
both short and long time scales. In III-V ferromagnetic semiconductors such as GaMnAs 
and InMnAs a small concentration of Mn ions is randomly substituted to cation sites so that 
the Mn-Mn spin coupling is mediated by the hole-ion p — d exchange interaction, allowing 
the generation of a ferromagnetic state with a Curie temperature of the order of 50 K . 
The magnetism can therefore be efficiently modified by controlling the hole density through 
doping or by excitation of electron-hole pairs with a laser pulse. In particular, unlike metals, 
in a regime of strong laser excitation total demagnetization can be achieved {4]. 

In the Zener model j^j], which was originally developed to describe the magnetism of tran- 
sition metals, the d shells of the Mn ions are treated as an ensemble of randomly distributed 
impurities with spin 5/2 surrounded by a hole gas or an electron gas. Unlike ferromagnetic 
metals, III-Mn-V ferromagnetic semiconductors offer the advantage of providing a clear dis- 
tinction between localized Mn impurities and itinerant valence-band hole spins, thus allowing 
the basic assumptions of the Zener theory to be satisfied. Based on this hypothesis, a few 
mean-field models have been successfully applied for modelling the ground-state properties 
of DMS nanostructures. In particular, within the framework of the spin-density-functional 
theory at finite temperature, relevant predictions of the Curie temperature have been ob- 
tained 6|j \]\. Ultrafast demagnetization in DMS is a phenomenon where the p — d exchange 
interaction causes a flow of spin polarization and energy from the Mn impurities to the holes, 
which is subsequently converted to orbital momentum and thermalized through spin-orbit 



and hole-hole interactions 



8|. Since energy and spin polarization transfer is a many-particle 



ying explanation of the ultrafast 
3- 



effect, the mean-field Zener approach cannot provide a satisf 

demagnetization regime that has been observed in DMS 

A phenomenological approach able to take into account this energy flux was given in 
[l(J] where a model based on three temperatures was derived. More recently, a study 

of the coupling of the electromagnetic laser field with the hole gas revealed the possibility 
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of an ultrafast demagnetization during the femtosecond optical excitation, due to light- 



hole entanglement 
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A model capable of describing the dynamics of carrier-ion s pin 



interactions is provided in 
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131 ] . This model generalizes the stationary theory of 10] 



and takes into account the picosecond demagnetization evolution which occurs in a strong 
excitation regime, but neglects the slow-in-time evolution of the spin dynamics. The mean 
p — d interaction is averaged out over the randomly distributed positions of the Mn ions. 

In this work we derive a dynamical model based on a many-particle expansion of the 
p — d exchange interaction based on the pseudo-fermion method. This formalism, originally 
developed by Abrikosov to deal with the Kondo problem [ijl], introduces unphysical states 
in the Hilbert space for which impurity sites are allowed to be multiply occupied. Following 
the work of Coleman [ijj], a suitable limit procedure is applied to our dynamical model in 
order to recover the correct physical description of the magnetic impurities. 

Our approach extends the Zener model beyond the usual mean-field approximation. It 
includes both the subpicosecond demagnetization dynamics and the slower cooling processes 
that restore the initial ferromagnetic order (which is achieved in a ns time scale). Moreover, 
in agreement with recent experimental results 9| our simulations show that, depending 
on the initial lattice temperature, a subsequent enhancement of the total magnetization is 
observed within a time scale of 100 ps. 



II. PSEUDO-FERMION FORMALISM 

We consider a volume V containing N h V holes with spin S = 1/2 strongly coupled 
by spin-spin interaction with N M V randomly distributed Mn impurities with spin S M = 
5/2. We assume that the exchange interaction between localized ions and heavy holes 
dominates over both the short-range antiferromagnetic d — d exchange interaction between 
the ions and the s — d exchange interaction between electrons in the conduction band and 
Mn ions (typical values for the s — d and the p — d interactions in a GaAs are 0.1 eV and 



1 eV respectively [161] V Furthermore, electron- hole radiative recombination, carrier-phonon 
interactions, and interactions leading to the hole spin-relaxation in the hole gas are included 
phenomenologically. The time evolution of the system is governed by the Hamiltonian 
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where a\ s (ojt, s ) is the creation (annihilation) operator of a hole with spin projection s 
and quasi-momentum k. In the parabolic band approximation the kinetic energy of the 
holes reads €k, s — E h — ^K- where E h is the valence band edge. The Kondo-like exchange 
interaction 7i p d is given by 

1 >-pd — y / J J m' ,m ' u s',s y u r],m' V< m k',s' k,s J e 

where the sum is extended over all indices, 7 is the p — d coupling constant, and <x, J are the 
spin matrices related to S h and S M respectively. The ion spin operator is represented in the 
pseudo-fermion formalism [14I. \i\ in which b\ jm (b v ^ m ) denotes the creation (annihilation) 
operator of a pseudo-fermion with spin projection m and spatial position R^. 

The 7i p d Hamiltonian reproduces the correct ion-hole exchange interaction provided that 
the ion sites are singly occupied, i.e., n v = X]m=-s M fyLm^m = ^ ^V- Following [u], 2] this 
constraint may be taken into account by adding a "fictitious" ionic chemical potential 

7Y A = ^ y Xrj fir) 

V 

to the original Hamiltonian and letting A,, go to infinity at the end of the calculation. 

The grand-canonical expectation value of a pseudo-fermion operator A related to the 
total Hamiltonian 7i + 7i x reads 

(A) x = ^Tr{ Pn e-^ x ^A} 

= IT E (K\pne~^ XvAv A\n™) , 

where Z x = Tr {p H e~ 13 ^ A """ }, p H = e~ m : , and (3 = l/k B T h , with k B the Boltzmann 
constant and T h the hole temperature. {^™} r = ^n\, . . . ,nf S +1 \ . . . ,ni 2S denotes 
all possible occupation numbers (= or 1) for r ion sites. Since each site has (2S M + 1) 
available pseudo-fermion states, the system will contain at most (2S + l)r pseudo-particles. 
The correct expectation value of the operator A is obtained using the limit X v — > 00 [l^] 

(A) = — lim d ZMA^ (1) 

where Z^ = \im {Zri} ^ d ^.% r and z n = e" /3A ". 

In the next section, we will show that the time-evolution of the spin of the ion-hole 
system may be expressed in terms of the expectation value of the pseudo-fermion operator 



fyr7,mV™»(l — b lim ,b v,m') with m 7^ m ' an d evaluated in the mean magnetic field S generated 
by the holes. We have the general relationship (which also applies when the system is driven 
far from equilibrium) 

l\ m (fymKmO- ~ b l,m>hm')) = A™ ( b l,mKm) A ■ ( 2 ) 
{A^j^oo \ / A {A,,}— >oo A 

When the system approaches thermal equilibrium, the quantity (b^^b^m) becomes the 
usual spin thermal distribution. Using Eq. (JT]) we obtain 

= Q , (3) 

-£■00 

where Z M = Q and Q = Q^^^ with Q = n^^ e 

In order to derive Eq. ([3]), we have used 

lim — Tr \p n e~^' VV n™| = Q V nfe~^ mn ™ = Q e~^ s m 

{^}->o dz\ ■ ■ ■ oz r 1 ' ' ^—f 1 

with p w = e" /37S ^. mmn '' ri and = b v>m . 

III. TIME EVOLUTION MODEL 

The Heisenberg equations of motion lead to a hierarchy of time evolution equations for 
the mean densities n h s = ^ ]Cfc(4, s a M) oo and n% = ^ X^(^,nAm)oo 

d[E t (4,, a ,,,)A] = (4) 

mi 



(5) 



with 



si ,m 



! 'ii m 'i 



In the last equation, the mean correlation function reads 

C , , = 1 7 r V,V,k[M iikx-k'JRr, (7] 

Lm > m i- S i> s i hVN h N M 2-^1 m'.mi^i.si ' ^ 

J7,fci,fc'i 
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where C^/^^g = (bL m >b Vtm ali s /Ofc, s )A- The time evolution equation of this quantity is given 
by 

^ ' = AEmf Cm',m,s',s (8) 

where the compact notations m = (r),m), s = (fc, s), £> = 6 m , 6 mi 6 m /6 m , H* = b^bmb^, b mi , 
A = a' a , a Sl al,a s , A 1 = ag,a s Og, a Sl have been employed. 

The mean-field contribution to the total energy is given by AEmf — 
7 [(s' -s)M + {m'- m) S] where M = N M Em=-s« m n ™ and S = Nh E£-5* s n s are 
the mean magnetic field generated by the ions and by the holes respectively. 

The use of Eq. (JS} combined with Eqs. (jl])-© leads to a non-Mar kovian time evolution 
of the macroscopic dynamical variables such as the density and the magnetization. By 
assuming an instantaneous spin-spin interaction, the Markov approximation can be easily 
recovered. For further details about the justification of the Markovian approximation in a 
DMS excited by a laser pulse we refer to [la ]. 



By using the Dirac identity [17[ e dt' = —ttM (e) — iKP\ where V denotes 

the principal value, the integration of Eq. (|Sj) with respect to the time leads to 

Cm',m,s',s = S ( A1E MF) 

mi,m'j,si,sj 

x 3 m[ , mi ■ <t s[iS1 (BA - B t A t ) x e #i-^i+C*'-^] , ( 9 ) 

where AEmf = £k' — £fc + AEmf- 

Since the matrix operators J • er are real, it is clear from Eq. ((6]) that the imaginary part 
gives no contribution to the equation of motion. 

The many-particle expansion of the correlation function C allows us to express Eq. Q in 
terms of the single-particle density matrix elements and n^f . By using the commutation 
rules of the creation and annihilation operators we obtain 



{BA— B t A t ) x = 5mi,m'£m ) m' 1 ^si,s'<$si 1 s 

((binbm -b ] m ,b m , J a\a s (l 
+ blb m ( 1 - bi,b m ,) ( 



al,a s / 



ala s — a[,a s ' 



A 
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Furthermore, as a closure hypothesis, we have assumed that the non-diagonal matrix el- 
ements of the density-like operators a\,a s and $ m ,b m with respect to the indexes 7] and k 
vanish. From the above approximations and using the definition (J3J) and Eq. we get 

Hi si,mj,mi 

5 (ae M f) (n* im , S)Sl - n^ mi:Sl)S J (10) 

where, 

k,kx,rj 

A similar expression can be found for ^2 W SljSlirriim in Eq. (JSJ). By using Eq. (J2J) we recover 
the fermionic limit of II, namely 

n^ limilS>sl = N M n%, (4,s*k,s (l - 4 1>Sl ^ llS1 ) ) oo • (12) 

In order to evaluate the time derivative of n(f, Eq. (1121) can be solved numerically. In the 
following paragraphs, we show that Eq. (fl2"l) may actually be further simplified. According 
to the Zener model the ground state of the system can be estimated by taking into account 
only the mean-field interaction between the holes and the magnetic ions. The hole gas 
experiences a mean magnetic field equal to M and in turn generates a mean field acting 
on the ions system equal to S. By converting the sum over k and k\ in Eq. (TT2"j) into the 
corresponding integral with respect to the energy variable E = we obtain 

TTOO ae mf r 

= N M n^,e~^ J f a h p(E)p(E - AE MF ) dE , (13) 

where /„ = (a^a s )oo (l - (4i a si)oo), h = ^^^^jj^i, and P denotes the hole density 
of states. 



In the limit 7S <C 7M <C we have 
TI 00 /9m 



V 2 



^j^3^e-^^(l-<). (14) 



In the next section we will validate this approximation by comparing the time evolution of 
the magnetization obtained by using either the approximate formula (j!4p or the exact one 
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(|T2|) . Finally, by inserting Eq. (|T4j) into Eq. (TTUT) we obtain 



5 M -1 



df 



2 ^ M R E - ™) s " + ™ + !) ( Z - A " 1/2 " ^Vi 2 ' 1/2 ) (15) 



m=-5 M 



dn 



M 



2£N h ( SM -om+1) (S M + am) (z~^ - Z^ 2 ) (16) 

cr=±l 

where 

AE MF 

Z S J = n%n h s (l-n h s ,)e ^ , (17) 



with f = 2vr7 2 ^v / 37T 2 A^ /l . 

IV. SPIN EVOLUTION IN DMS 

In order to study the time evolution of the mean magnetization of a GaMnAs/GaAs DMS 
heterostructure occurring after the interaction with a linearly polarized femtosecond laser 
pulse, we have applied our time dependent mode, constituted of Eos. & and flj» Based 
on the experiment of |9j, we consider a sample consisting of a 73 nm Gao.925Mn .o75As layer 
deposited on a GaAs buffer layer and a semi-insulating GaAs substrate. The background 
hole density is chosen to be 10 20 cm -3 . For the details of the chemical composition of the 
sample we refer to 9|. We assume that before the laser is turned on, the ion- hole system is at 
equilibrium with the phonon bath at the lattice temperature T L , so that the ground state can 
be well described by the Zener-type model described in The laser excitation generates a 
non-thermal electron-hole pairs distribution. By means of the Coulomb hole- hole interaction, 
the hole distribution undergoes a quasi- instantaneous thermalization (within a few tens of 
femtoseconds) towards a Fermi-Dirac distribution with temperature T h and a spin dependent 
chemical potential /i^ Q]- The increasing of the overall temperature T h of the hole gas is 
determined by assuming that the excess of energy of the hot photo-created particles (which 
is estimated as a fraction of the pump pulse energy) is redistributed among the total number 
of holes. The photo-created particles are approximately 2% of the background hole density 
9|. In particular, we consider an excitation by a monochromatic laser pulse tuned at the 
eue rgy E, aud having a pu m p fluence P, To estoate the energy E CI transferred rmtiaUy 
from the electromagnetic field to the kinetic energy of holes and electrons, following [12J , we 
assume that the fraction of the laser pulse energy imparted to the holes is 1/4 of the photon 
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FIG. 1: Temperature dependence of the normalized equilibrium ion magnetization for two hole 



densities. Here Nj 



5 x 10 19 cm 



(el + sl) 



energy. The total injected kinetic energy is thus E ex = r^E'f] with E[ = E\ — E g 
and el, el are the first eigenvalues of the valence and conduction bands. n^ x is the density 
of photo-created particles and rj is the ratio of kinetic energy absorbed by the electron gas 
which can be estimated within the spherical band approximation as rj = m^ H / (mjj + m^ H ) 
with m¥ H (mm) the effective mass of the heavy hole (electron) in the parallel direction 
of the sample. In Fig. [T]we show the normalized equilibrium ion magnetization M/iV M as 
a function of the temperature and parameterized by the hole density N h . In the figure, Nq 
indicates the initial density of holes, N£ x the excess of holes excited by the laser pulse, and 
N h = Nq+N^ x . After the laser excitation, the magnetic impurities strongly interact with the 
out-of-equilibrium hole gas by means of the 7i p d exchange interaction which redistributes the 
spin polarization from one system to another while conserving the total spin magnetization. 
Meanwhile, the itinerant hole spin is efficiently dissipated through spin-orbit interactions 
(tso ~ 100 fs) and relaxation of the total magnetization can be observed. Short-time spin 
relaxation of the holes is therefore an essential ingredient for explaining the observed time- 



12|. 



dependent changes of the magnetization in ferromagnetic semiconductors 

By means of a standard relaxation model, we include both the spin-orbit mechanism 
(or any other mechanisms leading to the hole-spin relaxation) and the other thermalization 
effects, such as the cooling of the kinetic energy of the excited holes driven by the phonons, 
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FIG. 2: Evolution of the normalized holes magnetization (S) and impurities magnetization (M), 
in the S - M plane, for T L = 70 K. 

and the radiative recombination of the electron-hole pairs. The corresponding equations 
read as follow 

^„ h „h „ h 

(18) 
(19) 
(20) 

where T h {t) and T L are the temperatures of the holes and the lattice, n^(n^ ,T h ) is the 
self-consistent quasi-static equilibrium hole spin distribution computed from the Zener-type 
model of which depends parametrically on the time- dependent ion magnetization M. 
The temperature relaxation rate t^ 1 



dn h s 






dt 


so 




dN h 


N h - N% 


dt 




dT h 




dt 





r p + r A p takes into account both the acoustic 



phonon scattering with tap = 200 ps and the optical phonon scattering with r OP = 1 ps 
for T h > 50 K and Top = oo for T h < 50 K [if]]. Eq. ffl9l takes into account the radiative 
recombination process characterized by a relaxation time trr = 400 ps [lj]]. 

In Fig. [2] we present the time evolution of the normalized magnetizations, where the 
vertical axis corresponds to S = S/N h and the horizontal axis to M = M/iV M . In agreement 
with the experiment of 9J, we consider a regime of small excitation (laser pump ffuence 
of 1 /iJcm -2 ) and a lattice temperature of 70 K. The point A represents the initial spin 
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FIG. 3: Time evolution of the total magnetization for different lattice temperatures: T L = 10 K 
(red line), T L = 30 K (green line), T L = 50 K (black line) and T L = 70 K (blue line). 

polarization which is suddenly shifted (instantaneously in our model) to point B. This is due 
to the laser excitation which abruptly enhances the hole density and consequently changes 
the normalization of S, so that S^O") = s ^ - and ^(O" 1 ") = s j° fe - . Our numerical simulations 
reveal the presence of different time evolution regimes: (i) < t < 50 fs: during this initial 
phase the magnetization evolution is nearly coherent (semi-coherent regime SC in Fig. [2). 
Indeed, since the photoexcited holes experience efficient spin-flip scattering with the localized 
Mn magnetic moments, a net spin polarization is transferred from the ion impurities to the 
holes leading to a significant increase of the hole spin polarization. Correspondingly, due to 
the large difference in densities between the two populations, only a small decrease of the ion 
magnetization is observed; (ii) 50 fs < t < 5 ps: the nonequilibrium hole spin polarization 
is efficiently dissipated via the spin-orbit coupling, which leads to a net decrease of the 
total spin magnetization (see also Fig. |3J). During this ultrafast demagnetization regime 
(UD) the kinetic temperature of the excited holes is still high; (iii) 5 ps < t < 350 ps: the 
hole distribution loses its energy via carrier-phonon scattering and the hole temperature 
decreases over the time scale tjj. When the Curie temperature is reached, the holes and 
ions spins begin to align, which allows the system to recover a ferromagnetic order. Since 
the total number of holes relaxes to its initial value N[} over a slower time scale trr ^> T£, 
a ferromagnetic state with an excess of holes can be reached, thus justifying a transient 
enhancement of the total magnetization (TE regime); (iv) 350 ps < t < 1.2 ns: finally 
the radiative recombination of the electron-hole pairs brings the system back to its initial 
configuration (RR regime). 

The time evolution of the total magnetization for different lattice temperatures is depicted 
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FIG. 4: Minimum (solid line) and maximum (dashed line) of the total magnetization for different 



lattice temperatures. 




100 200 300 400 500 600 
Time (ps) 



FIG. 5: Time evolution of the total magnetization at T L = 10 K. Full line: exact formula, Eq. 
(j 13|) : dashed line: approximated formula, Eq. (|14p . 



in Fig. [31 We see that the minimum of the total magnetization shifts to shorter times with 
increasing lattice temperature, in agreement with experimental findings. In Fig. H]we plot 
the excursion of the total magnetization for different lattice temperature: only for 45 K 
< T L < 78 K an enhancement of the total magnetization may be observed (4]. 

Finally, in order to validate the approximation of Eq. (1131) . we compare in Fig. [5] the 
time evolution of the total magnetization obtained by using either the approximate formula 
(I14p or by evaluating numerically the integral of Eq. (ITS]) . As can be clearly seen, a good 
agreement is obtained, justifying the use of the simplified expressions (fl~5|) and (TT6|) . 
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V. CONCLUSION 



In order to describe the strong spin-spin scattering regime observed in diluted magnetic 
semiconductors, we have derived a dynamical model that goes beyond the usual mean-field 
approximation. This model is based on the pseudo-fermion formalism and on a second-order 
many-particle expansion of the p — d exchange interaction, which is performed in terms of 
the single-particle density functions. At this level of description, this approach is similar to 
that of Ref. 12|, which was derived following a different perspective. Numerical simulations 



showed that our model is able to reproduce qualitatively - and to some extent quantitatively 
- the long-time evolution of the total magnetization after laser irradiation, as seen in recent 
experiments 9J. The early demagnetization observed in the experiments is explained as 
the result of a net flow of polarization from the ions to the holes, which is subsequently 
dissipated via spin-orbit coupling. Thus, the typical demagnetization time scale is mainly 
determined by the nonlinear coupling between the ions and holes spins, with a lower bound 
given by the spin-orbit time scale, tso ~ 100 fs. The demagnetization process cannot be 
faster than Tso, but can be significantly slower, depending on the lattice temperature. In 
addition - and in contrast to Ref. 12| - other slower processes (namely, holes thermalization 



and radiative recombination) were also included in the description, so that the global model 
encompasses time scales going from a few tens of femtoseconds to hundreds of picoseconds. 

We point out that the methodology developed in this work can be naturally extended 
to higher orders by using perturbative field-theoretic techniques. In particular, we plan to 
investigate third order dynamical processes, which were neglected here but may play an 
important role in the regime of higher photoexcitation energy. 

Finally, in the model used in this work only the heavy-hole band contribution to the 
exchange interaction was taken into account. The inclusion of a realistic band structure is 
currently under study. 
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